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Abstract 

The bound state energy eigenvalues and the corresponding eigenfunctions of the 
generalized Woods— Saxon potential are obtained in terms of the Jacobi polynomials. 
Nikiforov-Uvarov method is used in the calculations. It is shown that the results are 
in a good agreement with the ones obtained before. 
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1 Introduction 



Exact solution of Schrodinger equation for central potentials has generated much interes in 
recent years. So far, these potentials are the parabolic type potential pQ, the Eckart potential 
[JH2HS!; the Fermi-step potential [21 Ej, the Rosen- Morse potential [I], the Ginocchio barrier 
j3] , the Scarf barrier JGJ , the Morse potential [7j and a potential which interpolates between 
Morse and Eckart barriers fBJ. In addition, many authors have studied on exponential type 
potentials [01 EH HU El UHj and quasi exactly solvable quadratic potentials [TU UHl HE] ■ The 
exact solutions for these models have been obtained analytically. 

Recently, an alternative method known as the Nikiforov-Uvarov (NU) method has been 
introduced for solving the Schrodinger equation (SE). There have been several applications of 
SE with some well-known potentials fTJEHEI], Dirac, Klein-Gordon and Dumn-Kemmer- 
Petiau equations for a Coulomb type potential by using this method as well [2H [221 12H] ■ 

In the present work, the real- valued bound-state energies of the generalized Woods-Saxon 
potential are evaluated through the NU Method [21j by following the framework of quantum 
mechanics. This method is based on solving the time-independent Schrodinger equation by 
reduction to a generalized equation of hypergeometric type. Energy eigenvalues and the 
corresponding eigenfuctions are calculated exactly. The generalized Wood-Saxon potential 
as a shell model is selected. It can be used for describing metallic clusters in a successful way 
and for lighting the central part of the interaction neutron with one heavy nucleus [23| I26|. 

This paper is arranged as follows: In Sec. II we introduce the NU method. In Sec. Ill 
we apply the method to solve the Schrodinger equation with the generalized Woods-Saxon 
potential. Shapes of the potential and energy are studied in Sec. IV. In Sec. V we discuss 
the results. 

2 Nikiforov-Uvarov Method 

The NU method provides us an exact solution of non— relativist ic Schrodinger equation for 
certain kind of potentials |24j . The method is based on the solutions of general second order 
linear differential equation with special orthogonal functions [2Z] • For a given real or complex 
potential, the Schrodinger equation in one dimension is reduced to a generalized equation of 
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hypergeometric type with an appropriate s = s(x) coordinate transformation. Thus it can 
be written in the following form, 

r(s) + ^Voo + ^Mm = o (i) 

a a z {s) 

where a(s) and a (s) are polynomials, at most second— degree, and r (s) is a first— degree 
polynomial. To find a particular solution of Eq.Q by separation of variables, we use the 
following the transformation 

iP(s) = <f>(s)y(s) (2) 
This reduces Schrodinger equation, Eq.(l), to an equation of hypergeometric type, 

a(s)y" + r(s)y' + Xy = 0, (3) 

where <j>(s) satisfies 0(s)'/0(s) = ir(s)/a(s). y(s) is the hypergeometric type function whose 
polynomial solutions are given by Rodrigues relation 

^=w)^ l " n{SMS)l (4) 



where B n is a normalizing constant and the weight function p must satisfy the condition 

Op)' = tp- (5) 

The function 7r and the parameter A required for this method are defined as 



O- T 
7T = ± 



\ 



2 

— cr +ka (6) 



and 

A = k + ti'. (7) 

Here, ir(s) is a polynomial with the parameter s and the determination of k is the essential 
point in the calculation of 7r(s). Thus, in order to find the value of k, the expression under 
the square root must be square of a polynomial. Hence, a new eigenvalue equation for the 
Schrodinger equation becomes 

X = \ n = -nr' - n[n ~ l) a", (n = 0, 1, 2, ...) (8) 

where 

r(s) =t (s) + 2vr(s), (9) 
and it will have a negative derivative. 
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3 Generalized Woods— Saxon Potential 



The interactions between nuclei are commonly described by using a potential that consist 
of the Coulomb and the nuclear potentials. These potentials are usually taken to be of the 
Woods-Saxon form. As an example it can be given as the generalized Woods-Saxon potential 



V(r) 



C.e 



■r-Rp- 



1 + 



' r—Rg % 



2 ' 



(10) 



where Vo is the potential depth, Rq is the width of the potential and a is the surface thickness 
which is usually adjusted to the experimental values of ionization energies. In order to 
calculate the energy eigenvalues and the corresponding eigenfunctions, the potential function 
given by Eq. (JTUJ) is substituted into the radial part of Schrodinger equation: 

Vn Ce 2ar 



2m 



E + 



1 + qe 2ar (1 + qe 2ar ) 2 



if)(r) = 0. 



11) 



Here, some assignments are made in the radial Schrodinger equation such as R(r) = ip(r)/r, 
r — Rq = r and 1/a = 2a. In addition, q is deformation parameter and an arbitrary real 
constant within the potential. 

Now, in order to apply the NU— method, we rewrite Eq. (jllj) by using a new variable of 



the form s 



n 2ar 



d 2 i){s) 1 dtp(s) ^ 



m 



ds 2 



s ds 2h 2 a 2 s 2 



E + 



Vn 



Cs 



l-qs) (1-qsY 



ip{s) = 0. 



(12) 



By introducing the following dimensional parameters 



mE 



>0 (E<0), 



P 



mV 



{(3 > 0), 7 



mC 



' ~ 2h 2 a 2 v h " 2h 2 a 2 ^ " ' 2h 2 a 2 

which leads to a hypergeometric type equation defined in Eq.flTJ): 

1 



(7 > 0) (13) 



d 2 ip(s) 1 — qs dip(s) 



-eq 2 s 2 + (2eq - (3q - j)s + (3 - e ^(s) = 0. (14) 



ds 2 s(l — qs) ds s 2 (l — qs) 2 
After the comparison of Eq. (jl4}) with Eq.([T}, we obtain the corresponding polynomials as 



r(s) = l — qs, a(s) = s(l - qs), a (s) = -sq 2 s + {2eq — f3q - j)s + (3 - e. (15) 
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Substituting these polynomials into Eq.(jB|). we obtain 7r function as 

tt(s) = -y ± ^(g 2 + Aeq 2 - Akq) s 2 + A{(3q + 7 - 2eq + k) s + 4 (e - /?) (16) 

taking er'(s) = 1 — 2gs. The discriminant of the upper expression under the square root has 
to be zero. Hence, the expression becomes the square of a polynomial of first degree; 

(V + Aeq 2 - Akq) s 2 + A((3q + 7 - 2eq + k) s + 4 (e - 0) = 0, (17) 

and 

A = [4 + 7 - 2eg + A;)] 2 - 4 x 4 (e - /3) (q 2 + 4eg 2 - 4A;g) = 0. (18) 

When the required arrangements are done with respect to the constant k, its double roots 
are derived as k\ t2 = (f3q — 7) ± qJ (e — j3) (1 + y). 

Thus substituting, these values for each k into Eq.(fT5|) following possible solution is 
obtained for 7r(s) 



7T S 



2 2 



for A;=(/?g- 7 ) + g v /(e- / 9)(l + ^; 



(2v^5 + Jl + ^ )gs - 2v^^ 



for k=^q- 1 )-q^e-m + ^ 

(19) 



It is clearly seen that the energy eigenvalues are found with a comparison of Eq. (J7j) and 
Eq.fjHJ). From the four possible forms of the polynomial 7r(s) we select the one for which the 
function r(s) in Eq.® has a negative derivative. Therefore, the function r(s) satisfies these 
requirements, with 



t(s) = l-2qs-[ (2de -/3 + J1 + A 1 /q )qs -2de-/3 



t'(s) = -2q - (2 Je -/3 + dl + Aj/q )q. 



(20) 



Hence, the polynomial 7r(s) is computed from Eq.tjl6j) as 



gs 1 

7T S = 

v ; 2 2 



(2Je-P+ dl + A 1 /q ) qs -2de-f3 



From Eq.(jHJ) we also get 



A = (/3q - 7) - 9 " W + ~\-\ 



(2 v /^+ v /l + 4 7 /g )g 



(21) 



:22i 
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and also 



A = A r 



nq 



(2(1 + v /^3 ) + V / 1 + 47A? ) + n(n - l)q. 



(23) 



It is seen that the parameter e has the following form 

n 2 



■nq 



1 

16 



+ — + (l + 2n) 
Q 



P 2 



ft 



1 + % + (1 + 2n) 



^ + 2 



(24) 



Substituting the values of e and /? into Eq. (jl3|) and by using the transformation 2a = 1/a 
in Eq. ()22j) . one can immediately determine the energy eigenvalues E nq as 



E 



h 2 



nq 



2ma? I 16 



/ 8mCa 2 
1 + ^2— + (l + 2n) 



4 ^ m ^o" 2 



1 + 



8mCa 2 
7i 2 g 



+ (l + 2n) 



+ 



mV a 2 

IT 



(25) 



Here, the index n is non-negative integers with 00 > n > and the above equation indicates 
that we deal with a family of the generalized Woods— Saxon potential. Of course it is clear 
that by imposing appropriate changes in the parameters a and Vo, the index n describes the 
quantization for the bound states and the energy spectrum. In addition, if the parameter 
C in Eq. (|23j) is adjusted to zero, solution reduces to the form obtained for the standard 
Woods-Saxon potential without regarding the parameter q. 

Let us now find the corresponding wave eigenf unctions. Due to the NU— method, the 
polynomial solutions of the hypergeometric function y(s) depend on the determination of 
weight function p(s) which is satisfies the differential equation [a(s)p(s)]' = r(s)p(s). Thus, 
p(s) is calculated as 



qs) 



(26) 



where v 



1 + + 



■h 



Substituting into the Rodrigues relation given in Eq.(J3J), the 



eigenfunctions are obtained in the following form 



y riq (s) = B n (1 - qs) 



>-l) s ~2^y^f3 



d n 
ds n 



[1 - qs) 



(27) 



where B n is the normalization constant and its value is 1/n!. Choosing q — 1, the poly- 
nomial solutions of y n {s) are expressed in terms of Jacobi Polynomials, which is one of the 
orthogonal polynomials with weight function (1 — s) v ~ l s 2 ^ £ ~^ in the closed interval [0, 1], 
giving [constant] Pn v ^ e_/3 ' " ^(1 — 2s) 27J. By substituting tt(s) and a(s) into the expression 
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4>(s)'/4>(s) = tt(s)/<j(s) and solving the resulting differential equation, the other part of the 
wave function in Eq.(j2J) is found as 

M = (i - sTsV^e, (28) 

where \x = v/2 and again taking q — 1. Combining the Jacobi polynomials and <p(s) in 
Eq.pHJ). the s— wave functions are constructed as 

Ms) = A n sV^(l - s y^ +1 pi 2 ^ - 2s), (29) 

where A n is a new normalization constant. 

4 Potential and Energy Eigenvalue Shapes 

For the generalized Woods-Saxon potential in Eq. (jl(J|) . the potential shape can be determined 
by the parameters Vo, a and C. Therefore, we have adopted a generalized Woods— Saxon 
shape for the real part of the optical— model potential |.29j and investigated the scatter- 
ing phenomenon for this potential only. The empirical values found by Perey et al. [30J 
r = 1.285 fm, a = 0.65 fm and Vo ~ 40.5 + 13A MeV are used. Here, A is the mass 
number of target nucleus and R = tqA 1 ^ 3 . The shape of the generalized Woods— Saxon 
potential given by Eq.(10) is illustrated in Fig.l. 

Fig. 2 shows that the energy eigenvalues as a function of the discrete level n for different 
values of the parameter a. Of course it is clear that by imposing appropriate changes in 
the parameters a and Vq, the index n describes the quantization of the bound states and 
the energy spectrum. Some of the initial energy levels for q = 1 value are presented by 
choosing C — 10 MeV, A = 56 which is the geometric average of the target mass number 
44 < A < 72 [HI]. 

5 Conclusions 

The exact solutions of the radial Schrodinger equation for the generalized Woods-Saxon 
potential with the zero angular momentum are obtained. Eigenvalues and eigenf unctions 
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obtained from the real form of the potential are computed. So, the wave functions are 
physical and energy eigenvalues are in a good agreement with the results obtained by the 
other methods. Nikiforov-Uvarov method is used in the calculations. We have seen that 
there were some restriction on the potential parameters for the bound states within the 
framework of quantum mechanics. That is, the value of C from the potential parameters is 
increased for the constant value of the parameter a, it is seen that the depth of potential 
increases rapidly. Therefore, if all the parameters of potential remain purely real, it is clear 
that all bound energies E n with n > represent a negative energy spectrum. We also point 
out that the exact results obtained for the generalized Woods-Saxon potential may have 
some interesting applications in the study of different quantum mechanical and the nuclear 
scattering systems. 
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Figure Captions 



Figure 1: Variation of the generalized Woods-Saxon potential as a function of r. 

Figure 2: The variation of the energy eigenvalues with respect to the quantum number n 
with V 48 MeV. The curves correspond to the different values of the range parameter a. 
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Figure 1 








Figure 2 
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